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Abstract 

We study the regularity of the extremal solution of the semilinear biharmonic equation /3A^u — rAu = 
(^i^^yi on a ball B C R^, under Navier boundary conditions u — Au = on dB, where A > is a 
parameter, while r > 0, /3 > are fixed constants. It is known that there exists a A* such that for 
A > A* there is no solution while for A < A* there is a branch of minimal solutions. Our main result 
asserts that the extremal solution u* is regular (sup^ u* < 1) for N < S and /3, r > and it is singular 
(supg u* = 1) for N > 9, /3 > 0, and r > with ^ small. Our proof for the singularity of extremal 
solutions in dimensions A'^ > 9 is based on certain improved Hardy-Rellich inequalities. 

1 Introduction 

Consider the fourth order elliptic problem 

(3A^u - tAu = -p-^ in fl, 

0<u<l inn, (Ga) 

u = Au = on dft, 

where A > is a parameter, r > 0, /3 > are fixed constants, and C {N > 2) is a bounded smooth 
domain. This problem with /3 = models a simple electrostatic Micro-Electromechanical Systems (MEMS) 
device which has been recently studied by many authors. For instance, see |3], [5], [7], [5], [S], [TU], [H], [H], 
|13j . [I4] . and the references cited therein. 

Recently, Lin and Yang [16] derived the equation {G\) in the study of the charged plates in electrostatic 
actuators. They showed that there exists < A* < oo such that for A £ (0, A*) (Ga) has a minimal regular 
solutions u\ (sup^ ma < 1) while for A > A*, (Ga) does not have any regular solution. Moreover, the branch 
A — > u\(x) is increasing for each x £ B, and therefore the function u* = limA/A* u\ can be considered as a 
generalized solution that corresponds to the pull-in voltage A* . Now the important question is whether the 
extremal solution u* is regular or not. In a recent paper Guo and Wei '15] proved that the extremal solution 
u* is regular for dimensions A'^ < 4. In this paper we consider the problem (Ga) on the unit ball in M.^ : 

PA^u ^ tAu = -p-^ in B, 

< u < 1 in S, (Pa) 

u = Au = on dB, 

and show that the critical dimension for (Pa) is iV = 9. Indeed we prove that the extremal solution of (Pa) 
is regular (sup^ u* < 1) for N < 8 and /3, r > and it is singular (sup^ u* — 1) tor N > 9, (3 > 0, and t > 
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with ^ small. Our proof of regularity of the extremal solution in dimensions 5 < A'^ < 8 is heavily inspired 

by [4] and [6] . On the other hand we shall use certain improved Hardy-Rellich inequalities to prove that the 

extremal solution is singular in dimensions N > 9. Our improve Hardy-Rellich inequalities follow from the 

recent result of Ghoussoub-Moradifam [17] about Hardy and Hardy-Rellich inequalities. 

We now start by recalling some of the results from [T^ concerning {P\) that will be needed in the sequel. 

Define 

X*{B) sup{A > : (Pa) has a classical solution}. 
We now introduce the following notion of solution. 

Definition 1 We say that u is a weak solution of {G\), if < u < 1 a.e. in 57, (^i^l^yi £ L^( fl) and if 

n Jn (1 - uy 

Say that u is a weak super- solution (resp. weak sub-solution) of {G\), if the equality is replaced with > (resp. 
<) for (/3 > 0. 

We now introduce the notion of stability. First, we equip the function space Ti, := H'^{Vl) n Hq{^) = 
W^'^{n) n H^in) with the norm 



ll^ll = ( / [t\V^\' + m^A^]dx 



1/2 

Definition 2 We say that a weak solution u\ of{G\) is stable (respectively semi-stable) if the first eigenvalue 
Ml. of the problem 

2A 

- TAh + pA'^h - -h ^ nh mil, h^Ah^O on dn (1) 

(1 - u\y 

is positive (resp., nonnegative) . 

The operator /3A^u — tAu satisfies the following maximum principle which will be frequently used in the 
sequel. 

Lemma 1.1 ('151) Let u G L^iVl). Then m > a.e. in provided one of the following conditions hold: 

1. ue C^{Tl), [iA'^u - tAu >Q onVl, and u = An = on dfl. 

2. w(/3AV - rAip) dx>0 for all < ip e W^^^iQ) n H^iQ). 

3. u e W'^^'^(n), u ^ 0. An < on dB, and ^ [f3AuA(p + TWuV(p]dx > Q for all < Lp e W'^^'^{n) n 
Moreover, either u = or u > a.e. in fl. 

2 The pull-in voltage 

As in [B] and [1], we are led here to examine problem (Px) with non- homogeneous boundary conditions such 
as 

PA^u - tAu = in B, 

a <u<l in B, (Pa, a, 7) 

u = a, Au — J on dB, 

where a, 7 are given. Whenever we need to emphasis the parameters (3 and r we will refer to problem 
(-Pa,q.7) as {Px^p.T.a,^)- In this section and Section 3 we will obtain several results for the following general 
form of (Pa, a, 7) 

(3A'^u - tAu = in 57, 

a<u<l in ri, (GA,a,7) 

u = a, Au — J on dil, 
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which are analogous to the results obtained by Gui and Wei for {G\) in [TF. 
Let $ denote the unique solution of 

$ = a , A* = 7 ondn. ^ ' 

We will say that the pair (a, 7) is admissible if 7 < 0, a < 1, and sup^ $ < 1. We now introduce a notion 
of weak solution. 

Definition 3 We say that u is a weak solution of {P\,a,j), if a < u < 1 a.e. in D,, i^i]^y^y2 G L^{n) and if 

[{u- <^){(3A'^ - tA^) = a / e w^^^n) n H^{n), 

where $ is given in We say u is a weak super- solution (resp. weak sub-solution) of (Pa, 0.7); if the 

equality is replaced with > (resp. <) for (p >0. 

Definition 4 We say a weak solutionu of (Px^a.-y) is regular (resp. singular) if\\u\\oo < 1 (resp. \\u\\oo = Ij- 
We now define 



and 



A* (a, 7) := sup |a > : (Pa, a, 7) has a classical solutionj 
A* (a, 7) :— sup |a > : (Pa, 0,7) has a weak solutionj. 



Observe that by the Implicit Function Theorem, we can classically solve (Pa,q,7) for small A's. Therefore, 
A* (a, 7) and A* (a, 7) are well defined for any admissible pair (a, 7). To cut down on notations we won't 
always indicate a and 7. For example. A, and A* will denote the "weak and strong critical voltages" of 

(^A,a,7)- 

Now let U he a weak super-solution of (Pa, a, 7) and recall the following existence result. 
Theorem 2.1 fJ15f) For every < / G L^{il) there exists a unique < u G L^{fl) which satisfies 

u{j3 A^ Lp ~ T Aip) dx ~ I fipdx, 



for all if e W'^^'^iQ) n H^{n). 

We can introduce the following "weak" iterative scheme: uq = U and (inductively) let m„, n > 1, be the 
solution of 

(ii„-$)(/3AV-rA^) = A / ^ \fipeW^'\n)nH^{n) 



(1 - ur^-ir 

given by Theorem l2.1l Since is a sub-solution of (Pa, a, 7), inductively it is easily shown by Lemma [Ol that 
< Un+i < it„ < [/ for every n > 0. Since 

by Lebesgue Theorem the function u = lim u„ is a weak solution of (Pa, a, 7) so that a < u < U . We 
therefore have the following result. 

Lemma 2.2 Assume the existence of a weak super- solution U of (Px.a,-y) ■ Then there exists a weak solution 
u of (Pa, 0,7) so that a < u < U a.e. in fl. 

In particular, for every A £ (0, A*), we can find a weak solution of (PA,a,7)- In the same range of A's, this is 
still true for regular weak solutions as shown in the following lemma. 

Lemma 2.3 Let (a, 7) be an admissible pair and u be a weak solution of (Pa, a, 7)- Then, there exists a 
regular solution for every < < A. 
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Proof: Let e £ (0,1) be given and let u = (1 — e)u + e$, where $ is given in By Lemma ll.ll 

supQ $ < supjj u < 1. Hence 

sup u < (1 — e) + e sup $ < 1 , inf w > (1 — e)a + e inf <I> = a, 
and for every < (p e W^-'^{Cl) f) H^{^) there holds: 

(u-$)(/3AV-tA</7) = (1-e) / (m-$)(/3AV-tAv7) = (l-e)A ^ 
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(l-ft + e(<i>-l))2 - ^ ^ 7o 

Note that < (1 — e)(l — u) = I — u + e(<& — 1) < 1 — u. So u is a weak super-solution of (^'(i-£)3a,q,7) so 
that supu < 1. By Lemma 12.21 we get the existence of a weak solution w of (P(i_c)3;^ „ ^) so that a < w < u. 

In particular, supw < 1 and w is a regular weak solution. Since e G (0, 1) is arbitrarily chosen, the proof is 
n 

done. □ 

Lemma 12.31 implies the existence of a regular weak solution Ux for every A G (0,A*). Introduce now a 
"classical" iterative scheme: uq = and (inductively) m„ — Vn + ^, n > 1, where v„ G W*''^{n) fl (f^) is 
the solution of 

PA'^Vn — rAvn = PA'^Un — tAu„ = -TT in Q and Au„ =0 on dfl. (3) 

(1 - Un^iY 

Since w„ £ M^'*'^(ri) fl i?Q(r2), it„ is also a weak solution of ([3]), and by Lemma Fl. II we know that a < m„ < 

Un+i < C^A for every n > 0. Since supM„ < sup?7a < 1 for n > 0, we get that (1 — Un-i)^^ £ L^(r2) and 

n a 

the existence of w„ is guaranteed. Since w„ is easily seen to be uniformly bounded in iJ^(f2), we have that 
u\ := lim Un does hold pointwise and weakly in H^{fl). By Lebesgue theorem, we have that u\ is a 

n — >+oc 

radial weak solution of {P\) so that supuA < supC/A < 1. By elliptic regularity theory [1], u\ e C°°(f2) and 

o o 

wa = Aua = on dfl. So we can integrate by parts to get 



(3{A'^u\ — tAu\)(p dx — / u\{l3A'^ip — rAip) dx — X 



for every ip £ W'^''^{il.) fl i?g(ri). Hence, ma is a classical solution of (Pa) showing that A* = A*. 
Since the argument above shows that u\ < U for any other classical solution U of (P^j, a, 7) with /i > A, we 
have that ux is exactly the minimal solution and ux is strictly increasing as A j A*. In particular, we can 
define u* in the usual way: u*{x) = \wi ux{x). 

Lemma 2.4 A*(il) < +00. 

Proof: Let it be a classical solution of (Pa. a, 7) and let ("0, /^i) with Aip — on d^l denote the first cigcnpair 
of 13 A^ ~ tA in H^{n) n H^{n) with i; > 0.' Now let C be such that 

((ra - /37)9.V " /3a9,(A^)) ^ C [ ^. 
an Jn 

Multiplying (Pa, 0.7) by ^ and then integrating by parts one arrives at 

A 



n 



(1 -u)2 



Since ^ > there must exist a point x G $7 where _ ^,-\yi — iiiu{x) — C < 0. Since a < u{x) < 1, hence 



one can conclude that A < supq^„<]^(/xim + C)(1 — m)^, which shows that A* < +00. □ 
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In conclusion, we have shown the following description of the minimal branch. 
Theorem 2.5 A* G (0, +00) and the following holds: 

1. For each < A < A* there exists a regular and minimal solution u\ of {Px,a,'y)- 

2. For each x £ the map X t-* u\{x) is strictly increasing on (0, A*). 

3. For A > A* there are no weak solutions of {P\^a,-y)- 

3 Stability of the minimal solutions 

This section is devoted to the proof of the following stability result for minimal solutions. We shall need the 
following notion of weak solutions, which is an intermediate class between classical and weak solutions. 

Definition 5 We say that u is an Ti—weak solution of (Pa, a, 7) if u — ^ G H^{il) D i?Q(r2), < m < 1 a.e. 
in il, ^ G L^{^) and 



where $ is given by 0j. We say that u is an Ti—weak super- solution (resp. an Ti—weak sub-solution) of 
(^^,0.7) if for if > the equality is replaced with > (resp. <) and u > (resp. <), Au < (resp. >) on 



Theorem 3.1 Suppose that (0,7) is an admissible pair. 

1. The minimal solution u\ is stable, and is the unique semi-stable Ti~weak solution of {P\.a,'y)- 

2. The function u* := lim u\ is a well-defined semi-stable Ti—weak solution of (Pa*, 0,7)- 

3. u* is the unique Ti—weak solution of (-Pa*, a. 7), o,nd when u* is classical solution, then fj,i{u*) = 0. 

4. If V is a singular, semi-stable Ti—weak solution of (Pa, 0,7), then v — u* and A = A* 
The main tool is the following comparison lemma which is valid exactly in the class Ti. 

Lemma 3.2 Let (a, 7) be an admissible pair and u be a semi-stable Ti—weak solution of (Pa, 0.7)- Assume 
U is a Ti—weak super- solution of {P\^a,-y)- Then 

1. u < U a.e. in fl; 

2. If u is a classical solution and fJ.i{u) — then U — u. 

Proof: (i) Define w := u — U. Then by means of the Moreau decomposition for the biharmonic operator (see 
[l9j and there exist wi and W2 G H'^{fl) n Hq{D,), with w = wi -\- W2, wi > a.e., PA'^W2 — tAw2 < 
in the 7Y— weak sense and f3AwiAw2 + tVwi.Vu'2 = 0. Lemma 1 1 . 1 1 gives that W2 < a.e. in O. 





dn. 



Given < ^ G C^{n), we have 




where f{u) := (1 — u) 



^. Since u is semi-stable, one has 



A / f{u)wl< I /3(Awi)2 +t|Vu;iP = / /JAwAwi + rVw.Vwi < A / {f{u)-f{U)) 
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Since wi > w one has 

' f'{u)ww, < I {f{u)-f{U))Wi, 



I 

Jo. 



In 

which re-arranged gives 

' fwi > 0, 



where f{u) — f{u) — f{U) — f'{u){u ~ U). The strict convexity of / gives / < and / < whenever u . 
Since Wi > a.e. in fi, one sees that w <0 a.e. in fJ. The inequaUty u < U a.e. in Q is then estabhshed. 

(ii) Since m is a classical solution, it is easy to see that the infimum of Hi{u) is attained at some ip. The 
function is then the first eigenfunction of /3A^ — tA — ^j^^^^j in H^{n) n Hq{Q). Now we show that (p is 

of fixed sign. Using the above decomposition, one has (fi — pi + (p2 where ipi G H^{il) Ci Hq{Q) for i = 1,2, 
^1 > 0, f^(3AipiA(p>2 + T\/(p>i.V(p2 = and /3A^(/?2 — tA(P2 < in the H— weak sense. If p changes sign, 
then pi ^ and ip2 < in (recall that either 1^92 < or (^2 = a.e. in 17). We can write now 

^/3(A(<^i-^2))^+t|V(^1-^2)P-A/'(7/)(^1-<^2)^ 

= ^i(u) < 



JJ(pi-(p2y 

^/3(Ay.)^+r|V^|^-Af(^V 

< — -r 2 

Jn ^ 

in view of (pi(p2 < —<fi'P2 in a set of positive measure, leading to a contradiction. 

So we can assume > 0, and by Lemma [Ol we have 93 > in il. For < i < 1, define 

g{t)^ I l3A[tU+{l-t)u]Ap} + T\/[tU+{l-t)u].V(p-\ [ f{tU + {l-t)u)p>, 
Jn Jn 

where ip> is the above first eigenfunction. Since / is convex one sees that 

9{t) >X f [tfiU) + (1 - t)f{u) - f{tU + (1 - t)u)] ip>0 
Jn 

for every t > 0. Since 17(0) — and 

g'{0)= I (3A{U ~u)Alp + tV{U ~u).W^-\f'{u){U -u)p) = Q, 



we get that 

3"(0) = -A / f"{u){U^ufp>>Q. 



Since f"{u)(p > in fi, we finally get that U — u a.e. in il. □ 
A more general version of Lemma 13.21 is available in the following. 



Lemma 3.3 Let (0,7) be an admissible pair and 7' < 0. Let u be a semi-stable Ti—weak sub-solution of 
{Px.a,'y) with u = a' < a, Au — [3' > (3 on dfl. Assume that U is a 71— weak super- solution of (Pa, a, 7) with 
U = a, AU = (3 on dfl. Then U > u a.e. in il. 

Proof: Let u G H^{n) n Ho{fl) denote a weak solution of PA'^u - tAu = l3A^{u - U) - tA{u - U) in fl 
and u = Am = on d^l. Since u — u -\- U > and A{u — u -\- U) < on i9f2, by Lemma TLII one has that 
u > u — U a.e. in Vl. By means of the Moreau decomposition (see [12] and [5]) we write u as u ~ w + v, 
where w,v ^ i/p (il), w > a.e. in il, (3A'^v — tAv < in a W— weak sense and Jj-^ (3AwAv -f rVw.Vu = 0. 
Then for < (/3 e W^'^i^l) D H^{n), one has 



[ pAuAp + rVu.Vip <X [ (fin) - f{U))ip. 
Jn Jn 
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In particular, we have 

[ /3A£tAw + rVu.Vu; < A / {f{u)-f{U))i 
Jn Jn 

Since the semi-stabihty of u gives that 



A / f'{u)w'^< / /3(Au;)2 +r|Vwp = / /3AuAw + rVu.Vw, 
Jn Jn Jn 

we get that 

/ f'iu)w^ < I ifiu) - f{U))w. 
Jn Jn 

By Lemma 1 1.1 1 we have v <0 and then w > u > u — U a.e. in Q. So we obtain that 

0< / ifiu)- f{U)~f'{u){u~U))w. 
Jn 

The strict convexity of / imphes that U > u a.e. in fl. □ 

We need also some a-priori estimates along the minimal branch ux. 

Lemma 3.4 Let (a, 7) be an admissible pair. Then for every A e (0, A*), we have 



a {1-uxr - Jn (1-ua)2' 
where <i> is given by In particular, there is a constant C > independent of A so that 



for every A G (0, A*). 

Proof: Testing {P\,a^^) on ma - * e W^-'^{n) n Hl{n), we see that 

^ / Tr^2 - I (-1 -'^)\' + - '^)f)dx >2xf 

Jn (1 - uxY Jn Jn (1 " 

In the view of f3A^^ — rA$ = 0. In particular, for i5 > small we have that 

1 < 1 / ("A - ^ I f 1 



{\u^\>S} - Ux)^ J {\u^-'S>\>S} ~ ux)^ S^Jn{l-ux) 



by means of Young's inequality. Since for 5 small 

1 

-*|<A-} (1 - "a) 

for some C > 0, we get that 



J{W>. 



n 



[l-uxf 

f {t\Vux\'' + l3\Aux\^)dx = / (/3AuAA$ + rVuA.V$) + A / 
Jn Jn Jn ~ 



for some C > and for every A £ (0, A*). Since 

MA - * 



uxY 

< S f (Tl^uxf + muxf)dx + Cs + C ( f . ^ .3 
Jn \Jn (1 - ■ux)-^ 



in view of Young's and Holder's inequalities, estimate ^ is finally established. □ 

Proof of Theorem l3.lt (1) Since ||ua||oo < 1, the infimum defining fii{ux) is achieved at a first eigenfunc- 
tion for every A G (0, A*). Since A <—>■ ux{x) is increasing for every a; G 57, it is easily seen that A t-^ fJ.i{ux) is 
a decreasing and continuous function on (0, A*). Define 

A*, :== sup{0 < A < A* : ^i(ua) > 0}. 

We have that A*, = A*. Indeed, otherwise we would have /Lti(uA^^) = 0, and for every G ^A^^, A*), 
would be a classical super-solution of (-Pa,, ,0,7)- A contradiction arises since Lemma 13.21 implies = "a,,. 
Finally, Lemma [321 guarantees the uniqueness in the class of semi-stable 7Y— weak solutions. 

(2) It follows from ([4]) that u\ ^ w* in a pointwise sense and weakly in _ff^(f2), and j^j-r G L^{n). In 
particular, u* is a iJ^— weak solution of (PA*,a,7) which is also semi-stable as the limiting function of the 
semi-stable solutions {ua}- 

(3) Whenever ||u*||oo < 1, the function u* is a classical solution, and by the Implicit Function Theorem we 
have that = to prevent the continuation of the minimal branch beyond A*. By Lemma l3.2| u* is 
then the unique 7i— weak solution of (Pa*, 0,7)- 

(4) If A < A*, we get by uniqueness that v = ux- So w is not singular and a contradiction arises. Now, 
by Theorem I2.5r 3) we have that A = A*. Since w is a semi-stable H— weak solution of (Pa*, 0,7) and u* is 
a Ti.— weak super-solution of (Pa*, a, 7), we can apply Lemma 13.21 to get v < u* a.e. in fl. Since u* is also 
a semi-stable solution, we can reverse the roles of v and u* in Lemma [3.21 to see that v > u* a.e. in fl. So 
equality v — u* holds and the proof is done. □ 

4 Regularity of the extremal solutions in dimensions < 8 

In this section we shall show that the extremal solution is regular in small dimensions. Let us begin with 
the following lemma. 

Lemma 4.1 Let N >5 and (u*. A*) be the extremal pair of (Pa). If u* is singular, and he set 

T ■= {re {0,1) ■.us{r)>u*{r)} (5) 

is non-empty, where us{x) := 1 — C^lxj^ and Cs > 1 is a constant. Then there exists ri G (0,1) such that 
us{ri) > u*{ri) and Aus{ri) < Au*{ri). 

Proof. Assume by contradiction that for every r with us{ri) > u*{ri) one has Aus{ri) > Aw*(ri). Since F 
is non-empty and 

us{l)^l~C5 <0 = u*{l), 
there exists si G (0, 1) such that us{si) = u*{si). We claim that 

us{s) > u*{s), 

for < s < si. Assume that there exist S3 < S2 < si such that u*(s2) — us{s2), u*{s3) — 1/5(53) and 
us{s) > u*{s) for s G (s3,S2)- By our assumption Au^ > Au*{s) for s G {33,82) which contradicts the 
maximum principle and justifies the claim. Therefore us{s) > u*{s) for < s < si. Now set w := us — u* . 
Then w > on Bs^ and Aw < in i?si. Since w{0) — 0, by strong maximum principle we get w = on 
Bs^ ■ This is a contradiction and completes the proof. □ 

Theorem 4.2 Let N >5 and {u*,X*) be the extremal pair of {Px)- When u* is singular, then 

1-u* < C\x\^ in B, 

where C := (^)^ and A := . 
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Proof. For 6 > 0, define us{x) := 1 — C^lxl* with Cs '■— + (5)3 > 1. Since > 5, we have that 
Us G Hf^^(R^) and us is a ?i— weak solution of 



(3A'us-rAus = ^^^^ + ^rCs{N~h\x\-i in . 
We claim that us < u* in B, which will finish the proof by just letting (5^0. 

Assume by contradiction that the set F := {r e (0, 1) : us{r) > u*{r)} is non-empty. By Lemma |4TT] the set 

A := {r e (0, 1) : usir) > u*{r) and Aus{r) < Au*{r)} 
is non-empty. Let ri G A. Since 

us{l)^l-Cs <0^u*{l), 

we have that < ri < 1. Define 

a := u*(ri) < ^^(ri), 7 := Au*(ri) > Aus{ri). 

_ 4 

Setting us^n ~ ^ {us{rir) — 1) + 1, we see that us^n is a 7i— weak super-solution of {Px*+s\ p r^^r a' 7')' 
where ^ ^ 

a' := ''(a — 1) + 1, j' — r^j. 

Similarly, define u*^{r) = r^^u*{rir) - 1) + 1. Note that A^u* - aAu* > in B and An* = on 95. 
Hence by maximum principle we have Au* < in B and therefore 7' < 0. Also obviously a' < 1. So, 
(a', 7') is an admissible pair and by Theorem 13. ir 4') we get that (u*^, A*) coincides with the extremal pair 
of {P-^ p j,-2^ ^, in B. Also by Lemma [2?2] we get the existence of a week solution of {Px^+sx p r^^r a' 7')- 
Since X* + 6X> \* , we contradict the fact that A* is the extremal parameter of [P^ ^ ^.-'^^ a' j')- ^ 

Now we are ready to prove the following result. 

Theorem 4.3 //5 < iV < 8, then the extremal solution u* of {P)\ is regular. 

Proof. Assume that u* is singular. For e > define 1^5(2;) \x\~ and note that 

iV2(7V_4)2 



{A<ff = {Hn + C»(e))|x^^+2^ where Hn := 



16 



Given i] e C^{B), and since > 5, we can use the test function T]ip G Hq{B) into the stability inequality 
to obtain 

where 0(1) is a bounded function as e — > 0. By Theorem 14.21 we find 

„2 



and then 

2A / \x\-''+'' <{HN + 0{e)) [ |x|-^+2^ + 0(l) 
Jb Jb 

Computing the integrals on obtains 

2A < HN + 0{e). 

Letting e ^ we get 2A < H^. Graphing this relation we see that A^ > 9. □ 
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5 The extremal solution is singular in dimensions A^^ > 9 

In this section we will show that the extremal solution u* of (i^A,/3,T,o,o) in dimensions > 9 is singular for 
T > sufficiently small. To do this, first we shall show that the extremal solution of (-Pa. 1,0. 0,0) is singular 
in dimensions TV > 9. Again to cut down the notation we won't always indicate that f3 = 1 and r = 0. 
We have to distinguish between three different ranges for the dimension. For each range, we will need a 
suitable Hardy-Rellich type inequality that will be established in the appendix, by using the recent results 
of Ghoussoub-Moradifam [17]. 

• Case N > 16: To establish the singularity of u* for these dimensions we shall need the classical Hardy- 
Rellich inequality, which is vahd for all tp e H^{B) n Hq{B): 



T6 .L 



{AipY dx > ^— — / f-- dx. (6) 



Case 10 < < 16: For this case, we shall need the following inequality valid for all ip e H^{B) n Hq{B) 

^ (jV-2)^(jV-4)^ r ^ 

- 16 L(|,;|2__^|^|f + l)(|^|2_|^|f) 

(iV-l)(iV-4)2 f ^2 



• Case iV = 9: This case is the trickiest and will require the following inequality for all (p E {B)(^H}^{B) , 
which is valid for N > 7 



where where 



Au\^ > / W{\x\)u\ (8) 

B J B 



K{r) 



and 



(^(r) = r--+^ + 9r-'' + lOr - 20. 



The next lemma will be our main tool to guarantee that u* is singular for N > 9. The proof is based on an 
upper estimate by a singular stable sub-solution. 

Lemma 5.1 Suppose there exist A' > and a radial function u G H^{B) D Wl^^^{B \ {0}) such that 

A'u < 77^^ for < r < 1, (9) 
(1 — u)^ 

u(l) = 0, Aw|^=i==0, (10) 

u is singular, (11) 



2(3 



(l-u) 

for some /? > A'. Then u* is singular and 



I ,3 < / (A^)' for all ^EH\B)r\Hl{B), (12) 

Jb (1 - uY Jb 



X* < A' (13) 
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Proof. By Lemma [5751 we have (|13p . Let 4- < 7 < 1 and 



P 

« := C-^f'\ (14) 

and define u := 1 — a(l — u). We claim that 

u* <u m B. (15) 

To prove this, we shall show that for A < A* 

ux <u in B. (16) 

Indeed, we have 

,2.,^_„A2.-s - aA' a3\/ 



A^(u) =aA^(u) < 



(l-u)2 (l-u)2" 
By (fT3|) and the choice of a 

a^X' < A*. 

To prove it suffices to prove it for a^X' < X < A*. Fix such A and assume that (|15p is not true. Then 

A = {0<i?< 1 \uxiR) >u{R)}, 

in non-empty. There exists < i?i < 1, such that u\{Ri) > u*{Ri) and Au\{Ri) < Au*{Ri), since 
otherwise we can find < si < S2 < 1 so that ua(si) = u{si), u\{s2) = m(s2), u\{R) > u{R), and 
Aux{Ri) > Au*(i?i) which contradict the maximum principle. Now consider the following problem 

= 71 ^ ^ 

(1 — u)"^ 

u — u\{Ri) on dB 

Au — Aux on dB. 

Then ma is a solution to the above problem while m is a sub-solution to the same problem. Moreover u is 
stable since, 

A < A* 

and hence 

2A 2A* 2A' 23 

< ^TTZ 7^^—, < 



(1 - u)3 - a3(l - u)^ 7(1 - u)3 (1 - u)3 ■ 

We deduce m < ma in B^^ which is impossible, since u is singular while ux is smooth. This establishes (fT5 
From psp and the above two inequalities we have 

2A* 2A' 3 

< -7Z 7^ < 



(1-U*)3 - 7(l-1t)3 ^ (1-u) 

Thus 

inf (i?) -^^^ % >0. 

This is not possible if u* is a smooth solution. □ 
For any to > | define 

w„i := 1 - ajv,mr^ bN,mr'^, 

where 

TO(7V-hTO-2) _ 1(^-2/3) 

m{N + m-2)-l{N-2/3) ""'"^ r7i(iV + to - 2) - |(7V - 2/3) ' 

Now we are ready to prove the main result of this section. 
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Theorem 5.2 The following upper bounds on X* hold in large dimensions. 

1. If N > 31, then LemmaUl] holds with u W2, X'pj = 27A and l3 ^ > 27A. 

2. // 16 < iV < 30, then Lemmal5l\ holds with u W3, A'^ = - 1, /Sat = 

3. If 10 < N < 15, then Lemma \5.1\ holds with u W3, A'^ < (3^ given in Table\^ 
4- If N — 9, then Lemma \5.1\ holds with u := U'2.8> Ag :— 249 < [3g :~ 251. 

The extremal solution is therefore singular for dimensions N >9. 

Proof. 1) Assume first that N > 31, then it is easy to see that aN,2 < 3 and 2^ ^ 27A < We 
shall show that W2 is a singular 7i— weak sub-solution of {P)ajj ,a which is stable. Note that W2 G II^{B), 
G < ^2 < 1 in and 

A^W2<-^^^in B\{0}. 

(1 - 

So it;2 is a 7i— weak sub-solution of (-P)27a- Moreover, 

W2^l-\x\i + {aN,2 - - \xf) < 1 - \x\t 

Since 27A < we get that 

54A / < Fat / 7^^^^ < ff^ / rR ^ / 

for all ip E C^{B). Hence, W2 is stable. Thus it follows from Lemma [5TT] that u* is singular and A* < 27A. 
2) Assume 16 < < 30 and consider 

4 T 
W3 := 1 - ON.zr^ -1- OAT^sr . 

We show that it is a singular 7i— weak sub-solution of {Ph^ _ ) which is stable. Indeed, we clearly have 
< W3 < 1 a.e. in B, W3 G H'^{B) and G L^iB). Note that 



< sup- r^^;r^3 f ^ = Hn [ ^< I [A^r 

0<r<l (oat,™ - 6w,mr"~3)3 ' ./b 



Using maple one can verify that for 16 < A^ < 31 

A'w3<-^ ^ on (0,1). 

(1 - WsY 

Hence W3 is a sub-solution of (Ph^ By Lemma ISTD m* is singular and A* < — 1. 

3) Assume 10 < A' < 15. We shall show that W3 satisfies the assumptions of Lemma for each dimension 
10 < A^ < 15. Using maple, for each dimension 10 < A^ < 15, one can verify that inequality pT)) holds for 
A^ given by Table [1] Then, by using maple again, we show that there exists (3pj > A^ such that 

(A^-2)2(Ar-4)2 1 (A^- l)(Af-4)2 1 



16 {\x\^~^^\x\^+'){\x\^-\x\^) 4 \x\^{\x\^-\x\^) 

Wn 



> 



{l-W3f 
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Table 1- 






AT 
IN 


\ / 


Pn 


9 


249 


251 


iU 




OD / 


11 


405 


574 


12 


502 


851 


13 


610 


1211 


14 


730 


1668 


15 


860 


2235 


16 < iV < 30 


1 

2 


Hn 
2 


iV > 31 


27A 


Hn 
2 



The above inequality and improved Hardy-Rellich inequality ([31]) guarantee that the stability condition ((20l) 
holds for /Jat > A'. Hence by Lemma TS-ll the extremal solution is singular for 10 < < 15. The values of 
Xn and Pn are shown in Table [1] 
4) Let u:=W2.8- Using Maple on can see that 

A 2 249 . „ 
A^u < m B 

(1 — uy 

and 

502 

---^ < W{r) for all r e (0, 1), 

(1 — u[r))-^ 

where W is given by ([33|) . Since, 502 > 2 x 249, by Lemma [5.11 the extremal solution u* is singular in 
dimension iV = 9. □ 

Remark 5.3 It follows from the proof of Theorem \5.2\ that for N > 9 and ^ sufficiently small, there exists 
u e H^{B) n Wf^{B \ {0}) such that 

A^u - ^Au < TT^^ for < r < 1, (17) 
p (1 — u)"^ 

u(l) = 0, Am|,.=i = 0, (18) 

u is singular, (19) 

/ 7T^-^ ^ / (^^)' + ^l^^l' f'''^ ^ ^ ^'(^) ^ ^o(S), (20) 

JB U ^ J B P 



and 



where > > are constants. Indeed, for each dimension N > 9, it is enough to take u to be the 
sub-solution we constructed in the proof of Theorem \5.S\. j3'^ := (3^ , A' < A" < (3. If ^ is sufficiently small 

so that —^Au < ^iSj)-^ on (0, 1), then with an argument similar to that of Lemma \5.1\ we deduce that the 
extremal solution u* of (Px^p^rfifi) is singular. We believe that the extremal solution of {Px,p,T,o,o) is singular 
for all P,T > in dimensions N > 9. 



6 Appendix: Improved Hardy-Rellich Inequalities 

We now prove the improved Hardy-Rellich inequalities used in section 4. They rely on the results of 
Ghoussoub-Moradifam in [17j which provide necessary and sufficient conditions for such inequalities to hold. 
At the heart of this characterization is the following notion of a Bessel pair of functions. 
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Definition 6 Assume that B is a ball of radius R in R^, V,W & C^(0, 1), and -pT-^Ty-^dr = +00. Say 
that the couple (V, W) is a Bessel pair on (0, R) if the ordinary differential equation 

(Bv^w) y"{r) + + ^)y'{r) + ^y{r) = 

has a positive solution on the interval (0, i?). 

The needed inequalities will follow from the following two results. 

Theorem 6.1 (Ghoussoub-Moradifam [l7\ ) Let V and W be positive radial C^-functions on i?\{0}, 
where B is a ball centered at zero with radius R in (N > 1^ such that ^N-iy(^j.-^ d,r — +00 and 

r^~^V{r)dr < +00. The following statements are then equivalent: 

1. (y, W) is a Bessel pair on (0, R). 

2. ^gV{\x\)\Vip\'^dx> J^W{\x\)ip'^dx for all if e C^{B). 

Theorem 6.2 Let B be the unit ball in (N > 5). Then the inequality 

\Au\^dx > dx + {N-1) l^dx, (21) 

holds for allue C^{B). 

We shall need the following result to prove ([2T|) . 



Lemma 6.3 For every u S C^([0, 1]) the following inequality holds 



' |,.'(r)pr^-idr > /' f s^^, r^''dr - {N - 1)(«(1))^ (22) 

Proof. Let (p := r^"^+i — 2{n-i) ^^'^ ^(^) r^^^. Define ip{r) = u{r)/(p{r), r 6 [0, 1]. Then 

1 pi nl nl 

\u'{r)\'^k{r)dr = / \tl;{r)\^\cp'{r)\'^k(r)dr + / 2Lp{r)(p' (r)ij{r)i^' {r)k{r)dr + / \Lp(r)\^\tP' {r)\'^k{r)dr 
Jo Jo Jo 

1 .1 

|V(r)p(|(^'(r)|2fc(r) - {k^^'y{r))dr + / \^{r)\^\yj' {r)\^k{r)dr + 7^2(1)<^'(1)<^(1) 
^0 
1 



> / \^ir)\\\^'{r)\'kir)-ik^^'nr))dr + i;'il)ip\lMl) 
Jo 

Note that V^(l)^'(l)v5(l) = "^(1)^ = -{N - l)w^(l). Hence, we have 



'\u'ir)\'kir)dr > f -u\r/ + ^^^^"(0)^, „ „ i),2(,) (23) 
Jo ^ 

(24) 



Simplifying the above inequality we get ((22)) . □ 

The decomposition of a function into its spherical harmonics will be one of our tools to prove Theorem [621 
Let u e C^{B). By decomposing u into spherical harmonics we get 

u = S^gUfc where Uk = fki\x\)ipkix) 
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and (<Pfe(a;))fc are the orthonormal eigenfunctions of the Laplace-Beltrami operator with corresponding eigen- 
values Ck = k{N + k - 2), k > 0. The functions /?, belong to u e C°°([0, 1]), /^(l) = 0, and satisfy 
fkii") = 0{r^) and f'{r) — 0{r*'^^) as r ^ 0. In particular, 

ifia ^ 1 and /o - jy^^U-i /as. "^'^ = IvZ;^ i|x|=i <rx)ds. (25) 
We also have for any fc > 0, and any continuous real valued W on (0, 1), 

\Auk\^dx= [ {AM\x\)-Ck^^^fdx, (26) 



and 

' Wi\x\)\yuk\''dx= f Wi\x\)\yfk\^dx + Ck I W{\x\)\x\~^fldx. (27) 



Now we are ready to prove Theorem 16.21 We shall use the inequality 

/o W{r)\^r^-^dr > ^ ^^ r^^^dr for all x £ C\[0,1]) with x(l) = 0. (28) 

~ 2(]V-1) '"'^ 

Proof of Theorem I6.2t For all > 5 and fc > we have 



Nwn Jb A^wat 7s 







NcON 



{f^{r)fr^-^dr + {N-lf / (/^(r))\^-^dr 

Jo 

+4 f fl{r)r^-' + 2{N~l) f n{r)fl{r)r^~^ 

Jo JQ 

-2cfc f;:{r)h{r)r^''dr - 2ck{N - 1) mr)fk{r)r^~^dr. 
Jo Jo 

Integrate by parts and use (pS)) for fc = to get 

\Auk\^dx > (\raT)Yr^-Ur+{N -l + 2ck) [\fl{r)fr''-Hr (29) 
Jo Jo 

+ (2cfc(n - 4) + c^) r"-V.'(r)dr + {N - l){h{l)f 
Jo 

Now define gkif) = "^^^^ and note that guir) — 0{r''^^) for all fc > 1. We have 

'(/^W)V^-3 = f\g[ir)fr^-'dr+ f 2g,{r)g',{r)r^-^dr + fgl{r)r^~'dr 
Jo Jo Jo 

= [\g'k{r)rr^-'dr - (N - 3) f gUry'^-'dr 
Jo Jo 

Thus, 

> ^^^T^ /' , ^ r^-'dr - (iV - 3) /' /,^(r)r^-5dr (30) 
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Substituting ^Ck Jq {fl^{r)yr^ ^ in (1^ by its lower estimate in the last inequality ([50]) . and using Lemma 
we get 



J-^^N Jb 4 Jo ^ ~ 2(jV-^l) '' ^ ^ ~ 2(jV-l) '^ ^ 



(^-1) l\f'k{r)fr'''^dT + Ck{N-l) f\h{r)fr^''dr 
Jo Jo 







4 /n r2 - ^ T-f + 1 4 /n r2 - ^ rf + 1 

"'o r 2(Ar-i)'^ "'o r 2{N-i)^ 

+ (^-1) [\f'^{r)fT^-Hr + Ck{N-l) [\h{r)fr''-^dr 
Jo Jo 

The proof is complete in the view of (|27p. □ 
We shall now deduce the following corollary. 

Corollary 6.4 Let N > 5 and B he the unit hall in M^. Then the following improved Hardy -Rellich 
inequality holds for all Lp G H^{B) H Hq{B): 

(A^- l)(A^-4)^ r ip^ 
^ 4 Jb |a:|2(|:E|2 - |a;|f ) 

Proof. Let a := 2(n-i) ^^"^ ^(^) — note that 

K _ 2 a(Ar-2) ^ 2 

V r 2 1 — ar'^^-'^ ~ f 

The function y(r) = y^-f +2 — l is decreasing and is then a positive super-solution on (0, 1) for the ODE 

y +(^ + y)yM + -^2/ = o, 

where 



4(r2 — r 2 )(r2 — ar 2 +i) 
Hence, by Theorem 16 . 1 1 we deduce 



r+1 - ' 2 Wi3 (|a;|2 -a|x|f+i)(|a;|2 - ) 
for all e n i?oH-B). Similarly, for V{r) ^ ^ we have that 

Jb ~ 2 ' Jb |x|2(|x|2 - |x|f ) 

for all ip E H^{B) n Hg{B). Combining the above two inequalities with ([^H) we get ([51]) . □ 
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Corollary 6.5 Let N > 7 and B he the unit hall in M^. Then the following improved Hardy -Rellich 
inequality holds for all ip G H^(B) fl Hq(B): 

|Au|2 > f Wi\x\)u\ (32) 

B J B 

where 

.1^, ,33) 

K{r) = ^ ' ^ ' 

ip[r) 

and 

ip{r) = r"^+2 + 9r"2 + lOr - 20. 
Proof. Let a :— 2(n-i) ^'^'^ ^('') ■— — ■ Then ip is a sub-solution for the ODE 

r V V[r) 



where 



Hence by Theorem 16. II we have 



K(r) 



9 ■ 



Similarly 



where 



W ^ / (34) 



^ > / W,{\x\W. (35) 

B Fl 



K(r) 



Combining the above two inequalities with ((22|) we get improved Hardy-Rellich inequality (|32p. □ 

Acknowledgment: I would like to thank Professor Nassif Ghoussoub, my supervisor, for his valuable 
suggestions, constant support, and encouragement. I also thank C. Cowan for useful discussions. 



References 

[1] S. Agmon, A. Douglis, L. Nirenberg, Estimates near the boundary for solutions of elliptic partial differ- 
ential equations satisfying general boundary conditions. I, Comm. Pure Appl. Math. 12 (1959), 623-727. 

[2] T. Brau, A decomposition method with respect to dual cones and its application to higher order Sobolev 
spaces, preprint. 

[3] C. Cowan and N. Ghoussoub: Regularity of the extremal solution in a MEMS model with advection, 
preprint. 

[4] C. Cowan, P. Esposito, N. Ghoussoub, and A. Moradifam, The critical dimension for a fourth order 
elliptic problem with singular nonlinearity, submitted. 

[5] M. G. Crandall and P. H. Rabinowitz, Some continuation and variational methods for positive solutions 
of nonlinear elliptic eigenvalue problems, Arch. Ration. Mech. Anal., 58 (1975), pp. 207-218. 



17 



J. Davila, L. Dupaigne, I. Guerra, and M. Montenegro, Stable Solutions for the Bilaplacian with Expo- 
nential Nonlinearity, SIAM J. Math. Anal. 39 (2007) 565-592. 

P. Esposito, N. Ghoussoub, and Y. Quo, Compactness along the branch of semi-stable and unstable 
solutions for an elliptic problem with a singular nonlinearity, Comm. Pure Appl. Math., 60 (2008), pp. 
1731-1768. 

F. Gazzola and H.-Ch. Grunau, Critical dimensions and higher order Sobolev inequalities with remainder 
terms, NoDEA Nonlinear Differential Equations Appl., 8 (2001), pp. 35-44. 

N. Ghoussoub and Y. Guo, On the partial differential equations of electrostatic MEMS devices: Sta- 
tionary case, SIAM J. Math. Anal., 38 (2007), pp. 1423-1449. 

N. Ghoussoub and Y. Guo, On the partial differential equations of electrostatic MEMS devices II: 
Dynamic case, NoDEA Nonlinear Differential Equations Appl., 15 (2008), pp. 115145. 

Y. Guo, On the partial differential equations of electrostatic MEMS devices III: Refined touchdown 
behavior, J. Differential Equations, 244 (2008), pp. 2277-2309. 

Z. M. Guo and J. C. Wei, Hausdorff dimension of ruptures for solutions of a semilinear elliptic equation 
with singular nonlinearity, Manuscripta Math., 120 (2006), pp. 193-209. 

Z. M. Guo and J. C. Wei, Symmetry of nonnegative solutions of a semilinear elliptic equation with 
singular nonlinearity, Proc. Roy. Soc. Edinburgh Sect. A, 137 (2007), pp. 963-994. 

Z. M. Guo and J. C. Wei, Infinitely many turning points for an elliptic problem with a singular nonlin- 
earity, J. London Math. Soc, 78 (2008), pp. 21-35. 

Z. Gui, J. Wei, On a fourth order nonlinear elliptic equation with negative exponent, SIAM J. Math. 
Anal. 40 (2009), 2034-2054. 

F. Lin and Y. Yang, Nonlinear non-local elliptic equation modelling electrostatic actuation, Proc. R. 
Soc. A (2007) 463, 1323-1337. 

N. Ghoussoub, A. Moradifam, Bessel pairs and optimal Hardy and Hardy- Rellich inequalities, submitted. 

N. Ghoussoub, A. Moradifam, On the best possible remaining term in the Hardy inequality, Proc. Natl. 
Acad. Sci. USA 105 (2008) no 37, 13746-13751 

[19] J. -J. Morcau, Decomposition orthogonale dun espace hilbertien selon deux cones mutuellement polaires, 
C.R. Acad. Sci. Paris 255 (1962), 238-240. 



18 



